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A comparison is made of two different methods of calculating tran- 
sient t emp eratures for thick walls with arbitrary variation of heat- 
transfer coefficient and adiabatic-wall temperature . Although numerical 
calculations for special cases for which the exact solutions are avail- 
able show that both methods give satisfactory results. Hill's method 
(NACA Technical Note 4105) consistently gives nearly exact results with 
considerably less computing time, except for the case in which a tempera- 
ture profile through the thick skin is desired. For this case, 
Dusinberre's method (Trans. A.S.M.E., vol. 67 , no. 8) is much faster, 
though less accurate . 


INTRODUCTION 


.No exact analytical method is available for computing the transient 
temperature for the general case of thick walls . Various finite-difference 
methods have therefore been proposed and used to compute transient wall 
temperatures . With the increased importance of high temperature in air- 
craft structural design, an evaluation of the merits of representative 
methods is warranted. 

The existence of several other methods is fully acknowledged and 
no comprehensive comparison of all available methods is intended. It 
is intended merely to select two representative methods and to compare 
their results and computing times. The methods of Hill (ref. l) and 
Dusinberre (ref. 2) have been selected for this purpose. In this paper 
only the basic one-dimensional case has been considered. Hand calcula- 
tions with a desk computer were made rather than resorting to an elec- 
tronic computer, because in many engineering applications fast and direct 
answers are required, and in some problems, programing time on a computer 
becomes excessive. 

The first application of the finite-difference method for determining 
thick-wall temperatures is credited to Schmidt (ref. 5) • This method 
employs a ratio of incremental time to increment of distance into the 
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wall that is fixed "by the material properties of the wall. Dusiriberre 
has introduced an extension of Schmidt's method whereby the ratio of 
time increment to distance increment can be varied to introduce smaller 
time steps if desired. An increase of accuracy relative to Schmidt's 
method is thereby possible. Dusinberre's method includes Schmidt's 
method and can be reduced to it by the adjustment of a coefficient . 

Hill's method represents a considerably different approach to the 
thick-wall problem. Finite differences are taken only in the time var- 
iable, the equations used being already integrated with respect to 
distance . 

Following a description of both methods, sample problems are solved. 
Problems permitting exact solutions were chosen to make possible an 
evaluation of the accuracy of the methods . A time study of the methods 
was also made to determine the relative labor involved. An attempt is 
made to point out the areas of application wherein one method might be 
more advantageous than the other. 

It should be mentioned that Hill's method allows the outer- surf ace 
temperature to be readily determined if the time history of the inner- 
surface temperature is known. This fact makes Hill's method extremely 
advantageous to investigators in the field of aerothermodynamics where, 
generally, thermocouples are mounted on the inside surfaces of specimens 
and the heating rates and outer-surface temperatures are desired. The 
calculation of the outer- surface temperatures from the inner- surface 
temperatures can also be made with Dusinberre's method, but the process 
is much more laborious and less straightf orward . 


SYMBOLS 

c specific heat, Btu/(lb)(°F) 

F coefficient in Dusinberre's method 

G 
H 
h 
k 
K 

Z wall thickness, ft 


heat-capacity parameter, pcZ, Btu/(sq ft)(°F) 
heating-rate parameter, h5rt^/l6G 
heat-transfer coefficient, Btu/(hr) (sq. ft) (°F) 
diffusivity, K/cp, (sq. ft) /hr 
conductivity, (Btu) (ft)/(hr) (sq. ft)(°F) 
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M memory coefficient in Hill’s method 

m. step number 


P = cp(Ax) 2 / KB 
Q = h(Zflc) /K 


r. 

R 

t 

T 

5 

P 

e 

Ax 


radiation rate, Btu/(hr)(sq. ft) 
radiation-rate parameter, rBit 2 / l6G 
time from start of heating, hr 
temperature, °R 
time interval, hr 
weight density, Ib/eu ft 
memory coefficient 
distance increment, ft 


Subscripts : 

aw adiabatic wall 

i inner 

J distance increment number 

m step number 


OUTLINE OF PROBLEM AND METHODS 


The two methods were used to calculate the inner- and outer- surface 
temperatures of a thermally thick plane copper wall. Wall thicknesses of 
1/2 inch, 1 inch, and 3 inches were used. The walls were assumed to be 
insulated at the inner surface and their thermal properties were assumed 
to be constant. The given input function was a time history of adiabatic- 
wall temperature and heat- transfer coefficient. 
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Hill's Method 

Reference 1 gives a complete discussion and derivation of the equa- 
tions for Hill's method of computing transient temperatures of thick 
walls for any arbitrary variation of adiabatic-wall temperature and heat- 
transfer coefficient. The final equations, as presented in reference 1, 
are given here for convenience . 

Outer- surface temper ature . - The outer-surface temperature at the 
time mB is given by 


T™ = 


( HT aw)m + (®aw 


EE )m-l " ^ T m-1 “ M 3 T m-2 “ • * • “ *V T 1 “ ~ ^-1 

Mi + H m 


( 1 ) 


where 


For example, 


2 _ hSit 2 

(2) 

i6g 

R = rSjt2 

(3) 

16G 

G = pel 

oo 

(^aw)l + ( HT aw)o ~ R 1 “ R 0 

(5) 

M-i + H-j 


(®^aw)2 + (®aw “ *^)l ” ^2^1 " ®2 “ ®1 
Ml + Hg 


(6) 


Values of the memory coefficients M are taken from table I (most 
of which is reproduced from ref. l) . Interpolation in table I is avoided 
by working with a time increment 5 that results in a value of kB/Z 2 
listed In the table . 

Inner- surface temperature . - The inner-surface temperature at the 
time mB is given by 


T i,m = T m " ( 8 l T m + 6 2 T m-l + • • • + e m T l) ^ 

Values of 0 are taken from table I. 
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Dusiriberre's Method 


Dusinberre's method essentially consists of dividing the wall into 
a finite number of slabs and taking a heat balance for each slab . The 
relations used in calculating the trans lent- temperature time histories 
by this method are taken from reference 2 and repeated in this section. 


A remark to clarify the subscript notation Is in order. The first 
subscript on the temperature denotes the block for which the temperature 
is being calculated and the second subscript is the time at which the 
temperature is being calculated. On the averaging coefficients (the 
F- coefficients) the first subscript denotes the temperature used in the 
averaging process and the second subscript shows the block for which the 
temperature Is being calculated. Subscript i is the inner surface and 
subscript j is any intermediate block. The following sketch shows this 
notation: 


T 


aw ; m — 



*■ i 9 

A L, 

T l,m 

t 

T 2,m 
/L 

T j-l,m 

\ i 


T ,]+l,m 

a 

^i-l,m 

Ax — 

r 

"*• Ax 

If 

-Sr Ax-— 



m 


Outer-surface temperature .- For the outer surface. 


T l,m ~ F aw, l^aw, m-1 + F l,l T l,m~l + F 2,1^2,m-1 


( 8 ) 


where 


F n -^S = h~ — 
aw,l P k F 


(9) 


■p _ 2 

F 2,l = p 


( 10 ) 


F l,l = 1 ~ F aw,l " f 2,1 


Q, = 


h Ax 


(H) 


K 


( 12 ) 
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and P is some number such that 

P > 2 + 2Q (13) 

The parameter P can have any value which satisfies equation (l3), the 
size of P chosen dictating the time interval used, as shown by the 
relation 


g _ cp(Ax) 2 
KP 




Equations ( 13 ) and (l4) impose a maximum time increment for a given 
Ax but there is no limit to the minimum value of time increment which 
can be used by increasing the parameter P. Obviously, a time interval 
could be chosen and a value of P calculated from equation (l^) . If 
the P selected is large enough to be greater than 2 + 2Q for all 
values of heat-transfer coefficient, a constant S can be used which 
makes computing considerably faster and easier, the reason for this 
being that only F aw p and Fp p have to be computed for each time 
interval. The other F-coefficients are constant throughout the problem. 


Inner- surface tenrperatur e . - For the inner surface, 
•^i,m = F i-l,i^i-l,m-l * F i,i^l,m-1 


where the subscript i is the number of the last cube, and 


F — 2 

F i-l,i " p 


(15) 

( 16 ) 


F i,i " 1 " F i-1, i 


(IT) 


Intermediate wall temperatures .- The following equation permits a 
calculation of temperatures in the interior of' the wall: 

m _ T 3-l,m-l + (M - 2)T,3, m _p + Tj + i, m _p ,„ ON 

m C 8) 

Use of this equation allows temperature profiles through the slab 
to be calculated without any difficulty. Indeed, it is a consequence of 
this method that the entire temperature profile must be calculated, because 
the temperature at any point -in the slab depends on the temperatures on 
either side of it at the preceding time interval as well as its own tem- 
perature at the preceding time interval. This may be advantageous 
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"because Hill's method uses a series of cosine terms for the calculation 
of internal temperature profiles, which makes this calculation consider- 
ably more awkward. 


Numerical Calculations 

In the numerical calculations special caseB were chosen for which 
exact analytical answers could "be obtained, in order that the accuracy 
of each method might be determined. 

The transient- temperature histories of thick copper walls were com- 
puted by both methods. The following thermal properties were assumed: 

K = 227 (Btu) (ft)/(hr) (sq ft)(°F) 

c = 0.09192 Btu/(lb)(°F) 

p = 56 O lb/(cu ft) 

k = K/cp = I 4 -.I 4 I (sq. ft) /hr" 

Two different heating cases were considered. 

Case I .- The heat-transfer coefficient h was assumed to be 
100 Btu/ (hr) (sq ft)(°F) and to be held constant. The adiabatic-wall 
temperature was assumed to increase linearly with time, from a value 
of 0° F at zero time to a value of 10, 000° F at 10 seconds . The tran- 
sient temperatures of 1/2-, 1-, and 3- inch-thick copper walls were com- 
puted. The outer- and inner- surface temperatures for case I are given 
in figures 1 to 3- 

Case II .- The heat-transfer coefficient and adiabatic-wall temper- 
ature were assumed to have an arbitrary variation for the l/2-inch and 
the 3-Inch copper wall. Values of h are the same for both examples, 
whereas T aw - is slightly different. The values of h and T aw used 

in this case were obtained as follows : An outer- surface temperature was 

assumed and the heat input required to give this temperature was calcu- 
lated by an exact analytical method. This heat input was then used to 
determine values of h and T aw .. This was done so that a solution with 
h variable could be obtained from an analytical method. The tempera- 
tures calculated by the methods of Hill and Dusinberre were then compared 
with the original assumed temperatures . The values of h and T aw used 

are given in the following table: 
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Time, sec 

T °F 

x aw> r > 

for - 

h BtU 

1 / 2 - in. wall 

3 -in. wall 

7 (hr) (sq ft) (°F) 

0 

0 

0 

36.0 

1 

2,485 

2,484 

45.0 

2 

4,094 

4,088 

52.2 

3 

4,932 

4,915 

60.0 

4 

5,263 

5,227 

66 .6 

5 

5,387 

5,325 

69.0 

6 

5,356 

5,265 - 

66.6 

7 

5,107 

4,986 

60.0 

8 

4,335 

4,184 

52.2 

9 

2,769 

2,591 

45.0 

10 

297.5 

100 

36.0 


The results for case II are shown in figures 4 and 5- 

Computations . - For Dusinberre's method a choice must be made of the 
size of distance increments into which the wall is divided. One diffi- 
culty of the Schmidt or Dusinberre method for relatively thin walls is 
the small time step required by equation (l4), and therefore computing 
times which are long relative to the time required for thicker walls . 

To keep the computing times within reason, the l/2-inch wall was divided 
into only two increments . The increments selected for the different wall 
thicknesses are as follows: 


Wall thickness, in. 

Number of increments 

1/2 

2 

1 

2 

3 

6 


Hill's method uses the whole wall thickness, and thus removes the 
necessity of choosing an incremental thickness . 


RESULTS AND DISCUSSION 


In case I, wherein the adiabatic-wall temperature increases from 
0° F to 10,000° F in 10 seconds, typical results were as follows: 
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Case I 



Dusihberre 

Hill 

l, 

In. 


Maximum error 

Comp . 


Maximum error 

. 

Comp . 

8, sec 

°F 

Percent 
max . temp . 

time, 

min 

8, sec 

°F 

Percent 
max . temp . 

time, 

min 

1/2 

0.15 

11 

< 2 

183 

fi .4172 
\ .7086 

3 

1 

< 1 
< 1 

28 

50 

1 

{:© 

25 

18 

6 

4 

62 

103 

Jl.131 

l 

< 1 

26 

3 

.6 

15 

4 

98 

1.020 

2 

< 1 

32 


For the L/2-inch wall the error with Hill's method is about -1/4 
Anri the computing time about 1/5 that of Dusinberre's method. For the 
3-inch wall the error is about 1/8 and the time is about l/3 that for 
Dusinberre's method. In general, the maximum time interval was used in 
computing by Dusinberre's method. However, it has been found that a 
smaller 8 improves the accuracy, as shown in figure 2. Since the 
value of P used for Dusinberre's method was close to 2, the results 
for Schmidt's method would be similar. As shown by the results from 
alternate time increments, since Hill's method gives accurate results 
from substantially fewer steps, it is a waste of time to use very fine 
steps . For the thermally thinner walls - for example the l/2-inch wall 
the small 8 required to satisfy the Dusihberre (or Schmidt) relations 
necessitates a long computation time . 

In case IT, wherein the adiabatic-wall temperature both rose and 
fell in a 10-second period while the heat-transfer coefficient varied 
also, the results were: r 


Case II 



Dusihberre 

Hill 

l, 

in. 


Maximum error 



Maximum error 

Computing 
time, rain 

8, sec 

°F 

Percent 
max . temp . 

6, sec 

°F 


1/2 

0.15 

2 

< 1 

128 

0.71 

0-5 

<1 

60 

3 

.6 

3 

2 

57 

1.020 

• 5 

< 1 

34 
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In this case the maximum error in the Dusihberre method was only 
2° F for the 1/2- inch wall. However, for an equal error, Hill's method 
required only l/5 as long. For the 3 -inch wall the times were about 
twice as large for Dusihberre' s as for Hill's method. On the inner sur- 
face the error in Dusihberre 's method, while small, was a substantial 
percentage of the rise of the inner-surface temperature. 

Dusihberre ' s method is well suited to obtaining the temperature 
distribution through the wall since in all cases the distribution is 
obtained as a necessary consequence of the computation. Such distri- 
butions are shown for case I for the 3- inch- thick wall in figure 6 . At 

6 seconds a temperature distribution computed from an exact formula from 
reference 1 is shown for comparison. The maximum error is 12° F or about 

7 percent. It is also possible to compute the temperature distribution 
by Hill's method. This was done by the equation of appendix C of refer- 
ence 1. The results are in almost perfect agreement with the exact theory. 


CONCLUDING REMARKS 


From the examples presented herein, as well as from other examples 
presented in reference 1 , it appears that for any reasonable step size 
Hill's method is, practically, an exact method. On the other hand, 
Dusihberre' s method, with reasonable step sizes, gives a good approxi- 
mation. The same statement applies to the Schmidt method. Hill's method 
is also substantially faster than Dusihberre 's method (or Schmidt's 
method) if only the two surface temperatures are required. If tempera- 
ture distributions through the wall are required. Hill's method is slower 
but practically exact . Either method is suitable for machine calcula- 
tions. Exact (classical) methods are not available except for special 
cases . 


Langley Aeronautical Laboratory, 

National Advisory Co mmi ttee for Aeronautics, 
Langley Field, Va., May 1 6 , 1958 . 
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TABLE I.- VALUES QF MEMORY CQEFETCIHSTS 




(a) Values of M 


\ . ~ r. O 

\ko/ L*- 

0.01 

0.02 

O.05 

0.1 

0.2 

0.5 

1,0 

2.0 

5.0 

1 

0.09281491 

0.13125147 

O.20752122 

0.29347746 

0.41495581 

0.64728217 

0.8568372 9 

1.02954122 

-1.15145422 

2 

-.01594178 

- .02252780 

-.03560872 

- .05052075 

-.O78366OO 

-.22817193 

-.51142943 

-.82682879 

-1.06920825 

5 

- .02682329 

-.03794254 

-.06002358 

- .08722726 

-.14861514 

-.30052435 

-.31662004 

-.20127594 

-.08224562 

4 

-.00957530 

-.01354153 

-.02166582 

-.03696546 

-.07338883 

- .08405219 

-.02634648 

-.00142616 

-.00000056 


- .00 553400 

_ 0078968! 

_ m *onnofl 

_ noA oRnArr 

Ail llg*gl;l; 

.. 6ViJi li T^rOi 


/wwi rv^C. 

TVWVWkA 

6 

-.00373390 

- .00528488 

-.00983262 

-.02025561 

-.02724662 

-.00712803 

-.00018948 

-.00000007 


7 

-.00274116 

-.00389295 

-.00809228 

-.01579610 

-.01663398 

-.00207577 

-.00001607 



8 

-.00212389 

-.00504410 

-.00695594 

- .012338&7 

-.01015501 

-.00060449 

-.00000136 



9 

-.00170888 

-.00249269 

-.00608368 

-.00964054 

-.00619962 

-.00017604 

-.00000012 



10 

- .00141415 

-.00211866 

-.00535620 

-.00753254 

-.00378485 

-.00005126 

-.00000001 



11 

-.00119613 

- .OOI8553I 

-.00472749 

-.00588550 

- .00231065 

- .00001493 




12 

- .00103005 

-.00166301 

-.00417648 

-.00459860 

-.00141064 

-.00000435 




13 

-.00090099 

-.00151704 

-.00369098 

- .00359308 

-.00086120 

-.00000127 




14 

-.00079666 

- .00140173 

- .00326234. 

-.00280743 

-.00052576 

-.00000037 




15 

-.OOO71651 

-.00150718 

-.00288362 

-.00219357 

-.00032097 

-.00000011 




16 

-.00064769 

-.00122688 

-.00254890 

-.00171393 

-.00019595 

-.00000003 




17 

- .00059289 

-.00115664 

- .00225306 

-.00133917 

-.00011963 

-.00000001 




18 

-.00054714 

- .00100^86 


ooio46 vs 

-.oonrrzvra 





19 

- .00050882 

-.OOIO3656 

-.00176041 

-.OOO81756 

-.00004459 





20 

- .00047652 

-.00098373 

- .00155610 

-.00063800 

-.00002722 





21 

-.00044898 

-.00093449 

-.00137549 

-.00049912 

-.00001662 





22 

-.00042534 

-.OOO88826 

-.00121584 

-.00038998 

-.00001015 





23 

-.00040474 

-.00084473 

-.00107473 

-.00030471 

-.00000619 





24 

-.00038678 

-.00080357 

-.00094999 

-.00023809 

-.00000378 





25 

-.00037081 

-.00076453 

-.00083973 

-.00018602 

-.00000231 





26 

-.00035655 

-.OOO72753 

-.00074227 

- .00014535 

-.00000141 





27 

-.00034370 

-.00069238 

-.00065612 

-.00011357 

-.00000086 





28 

-.00033195 

-.00065894 

-.00057997 

-.00008874 

-.00000052 





29 

-.00032118 

-.00062716 

-.00051265 

-.00006933 

-.00000032 





30 

-.00031123 

-.00059693 

-.00045315 

-.00005417 

-.00000020 






f 


NI VOW 



(b) Values of 6 


1 

0.01 

0.02 

0.05 

0.1 

0.2 

0.5 

1.0 

2.0 

H 

1 

1.00000000 

O.99999990 

0.99956261 

0.98873107 

0.92596579 

0.69945338 

0.45623848 

0.24814437 

0.09999955 

2 

-.00000019 

- .00019242 

-.02166158 

-.12553131 

-.31355030 

-.48642965 

- .41618822 

-.24630212 

-.09999908 

5 

- .00002495 

- .00333067 

-.06963678 

-.17636952 

-.23787104 

-.15098863 

- .03665380 

-.00182900 

-.00000045 

4 

-.00053480 

- .OH95026 

-.09012900 

-.14882949 

-.14587796 

- .04396967 

- .00310842 

-.00001515 

-.00000000 

5 

- .00144190 

-.02167827 

-.08988005 

-.11749067 

-.O8906596 

-.01280454 

-.00026361 

-.00000009 


6 

-.00544274 

-.02904777 

-.08284990 

-.09193126 

-.05437468 

-.00372885 

-.00002236 

-.00000000 


7 

-.00595604 

-.OO854572 

- .03361672 

-.07435643 

-. 071644 Q 9 

-.03319564 

-.00108589 

-.00000190 



8 

-.03599626 

-.06609621 

- .056136W 

-.O2O265&7 

- .00031622 

-.00000016 


r 

9 

- .01092556 

- .03687745 

-.05854671 

- .04386205 

-.03237227 

-.00009209 

-.00000001 



10 

u 

-.01296009 

-.01461454 

- .05679326 
-.03611078 


-.03427135 
- .02677771 

-.00735325 
- .00461124 

-.00002682 

-.00000781 

-.00000000 



12 

- .01590635 

- .03507182 

- .04048318 

-.02092260 

-.00281516 

-.00000227 




13 

- .OI68767I 

- .03383231 

-.03578599 

- .01634774 

-.OOI71865 

-.00000066 




111- 

-.01757368 
- .01004424 

-.03249148 

-.03163302 

- .01277320 

- .00104923 

-.00000019 




15 

-.03im88 

- .02796174 

-.00998026 

-.00064055 

-.00000006 




16 

- .01833036 

-.02973243 

-.02471644 

- .00779802 

-.00039106 

-.00000002 




17 

-.01846822 

- .02837688 

-.02184780 

-.00609293 

-.00023874 

-.00000000 




18 

-.01848810 

- .02705950 

-.01931206 

-.00476067 

- .00014575 





19 

- .01841482 

-.02578770 

-.02456606 

-.01707064 

- .00371972 

-.OOOO8898 





20 

-.01826876 

-.01508936 

-.00290638 

-.00005432 





21 

-.01806620 

-.02339609 

- .01333804 

-.00227088 

-.00003316 





22 

-.01782041 

- .02227782 

-.01178998 

-.00177434 

-.00002025 





23 

-.01754197 

- .02121042 

- .01042159 

-.00138637 

-.00001236 





24 

-.01723931 

-.02019254 

-.00921203 

-.OOIO8323 

- .OOOOO755 





25 

- .01691916 

-.01922244 

- .00814285 

-.00084638 

-.00000460 





26 

-.OI658696 

- . 01829^9 

-.00719776 

-.00066131 

-.00000281 





27 

-.01624686 

-.01741809 

-.00636236 

-.00051671 

-.00000172 





28 

-.01590228 

- .01658000 

-.00562392 

-.00040373 

-.00000105 





29 

- .01555588 

-.01578204 

-.00497119 

-.00031545 

-.00000064 





30 

-.01520973 

-.01502235 

-.00439422 

-.00024648 

-.00000039 
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Figure 1.- Case I. Temperatures of l/2-inch. copper wall surfaces. 
Adiabatic-wall temperature varies linearly from. 0° to 10,000° in 
10 seconds; h. = 100 Btu/(hr)(sq. ft)(°F). 
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Figure 2.- Case I. Temperatures of 1-Inch. copper wall surfaces. 
Adiabatic-wall temperature varies linearly from 0° to 10,000° in 
10 seconds; h = 100 Btu/(hr)(sq. ft)(°F). 
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Figure 3.- Case I. Temperatures of 3-inch copper wall surfaces. 
Adiatatic-wall temperature varies linearly from 0° to 10,000° in 
10 seconds; h = 100 Btu/(hr) (sq. ft)(°F). 
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Figure ll-.- Case II. Temperatures of l/2iinch copper wall heated 
according to assigned history of h. and T aT , r . 
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Figure 5.- Case II. Temperatures of 3-inch copper, wall heated according 
to assigned history of h and T aw . 
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Figure 6.- Case I. Temperature profiles through 5-inch copper wall. 
Adiahatic-wall temperature varies linearly from. 0° to 10,000° in 
10 seconds; h = 100 Btu/(hr)(sq. ft)(°F) . 
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